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We study the pion form factor in a broad range of spacelike momentum transfers within the 
local-duality version of QCD sum rules. We make use of the recently calculated two-loop double 
spectral density of the (AVA) correlator including O(l) and 0(ct 3 ) terms, which allows us to give 
predictions for the pion form factor and to study the interplay between the nonperturbative and 
perturbative contributions to the pion form factor without any reference to the pion distribution 
amplitude. Our results demonstrate the dominance of the nonperturbative contribution to the form 
factor up to relatively large values of the momentum transfer: namely, the nonperturbative O(l) 
term, which provides the 1 / Q 4 power correction, gives more than half of the pion form factor in the 
region Q 2 < 20 GeV 2 . 



PACS numbers: 11.55.Hx, 12.38.Lg, 13.40.Gp 



INTRODUCTION 



The study of the interplay between perturbative and nonperturbative physics in exclusive processes, and, in par- 
ticular, in the pion form factor which we discuss in this paper, has a long history. At asymptotically large Q 2 , the 
pQCD factorization formula [l[ gives the pion form factor in terms of the scale-dependent pion distribution amplitude 
(DA) of leading twist w (m, Q 2 ): 



8 TOs (Q 2 )/ 2 
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The DA, as obtained from the pQCD evolution equation, has the form 



6m(1 



Respectively, at asymptotically large Q 2 a direct prediction of pQCD reads pi Q 

Q 2 F 7T (Q 2 ) = 8ira s (Q 2 )f 2 . 



(1.1) 
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Sublcading logarithmic and power corrections to this formula should be taken into account at large but finite Q 2 . 
This is, however, a very difficult task. There are two competitive scenarios for the pion form factor at intermediate 
momentum transfers: 

The first scenario (A) is based on the assumption that power corrections are negligible in the region Q 2 > 3-5 
GeV 2 . The form factor is then given by the pQCD factorization formula (|1.1|) with the pion distribution amplitude at 
low normalization scale, which in this scenario turns out to have a double-humped "camel" shape with an enhanced 
end-point region very different from its asymptotic form (|1.2[) . This scenario, complemented by the analysis of 
Sudakov double logarithms provides the basis for the perturbative QCD approach to form factors at intermediate 
momentum transfers. 

In the second scenario (B), which we consider to be more realistic, the form factor is dominated by the nonper- 
turbative contributions up to rather high values of Q 2 , with the perturbative contribution remaining relatively small 
[H, 0] . The end-point behaviour of the DA at a low normalization scale is then similar to that of the asymptotic DA 
(jl.2|) . This scenario is supported by the fact that the soft contribution to the form factor alone can reproduce the 
pion form factor to a good accuracy for Q 2 up to several GeV 2 @, 0, H, OH, [HI ■ 

In [12j |. making use of the constituent quark picture, and in [13j, within light-cone sum rules, the pion form factor 
was analyzed by taking into account the nonperturbative O(l) contribution and the radiative 0(a s ) corrections. The 
form factor at intermediate Q 2 turned out to be sensitive to the details of the pion wave function — the Bethe-Salpeter 
wave function in pjj and the pion light-cone DA in [l3j] . Scenario B was favoured by these results. 

Unfortunately, the data on the pion form factor for Q 2 > 2 GeV 2 are not sufficiently precise, leaving room for 
speculations about the details of the pion DA at low normalization scale and, respectively, on the relative weights of 
the soft and the hard contributions to the pion form factor. 
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Therefore, it seems interesting to address the problem without a direct reference to the pion DA. The local-duality 
version of three-point QCD sum rules Q provides this opportunity. 

The local-duality sum rule is the Borel sum rule in the limit of an infinitely large Borel parameter. For the relevant 
choice of the pion interpolating current, the condensate contributions to the correlator vanish in this limit and the 
pion observables are given by dispersion integrals via the spectral densities of purely perturbative QCD diagrams. 
The integration region in the dispersion integrals is restricted to the pion "duality interval" . 

In this Letter we apply a local-duality sum rule to the pion form factor, making use of the recently calculated two- 
loop double spectral density of the pion form factor for massless quarks [14l . |15| . Such an approach has the following 
attractive features: (i) it is applicable in a broad range of momentum transfers starting from low to asymptotically 
large values, and (ii) it does not refer to the pion distribution amplitude. Therefore, it allows us to study in a relatively 
model-independent way the interplay between perturbative and nonperturbative dynamics in the pion from factor. 



2. SUM RULE 



We shall consider the pion form factor in the chiral limit of massless quarks and a massless pion. Let us recall 
well-known results for Borel sum rules: The sum rule for the pion decay constant is obtained from the OPE for the 
two-point function and reads [l6| 



-so 
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where p(s) = -r- (l + ^) + 0(a 2 ) is the perturbative spectral density. 

The Borel sum rule for the pion form factor is obtained from the OPE for the three-point function and reads @, 0] 



f 2 F„(Q 2 ) =T(Q 2 ,M 2 ,M 2 \s ) 



(a s G 2 ) 4ira s (qq) 2 
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Here, T(Q 2 , M 2 , M 2 \sq) is the perturbative contribution, which is obtained by the following procedure: One calculates 
the double Borel transform of the (AVA) correlator 
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T(Q 2 ,M 2 ,M 2 ) = — / d Sl d S2 exp(- Sl /M 1 2 )exp(- S2 /M|) A^(Q 2 , Sl , s 2 ) + a s A«(Q 2 , s 1; s 2 ) 



(2.3) 



and restricts the integration in the si-s 2 plane to the pion duality region. One then sets M 2 ~ * 2M 2 , Mf — > 2M 2 
and compares the two- and three-point sum rules for the same values of the Borel parameter M 2 . 1 The function 
A(°)((3 2 , si, s 2 ) is well-known, whereas A^ 1 ) (Q 2 \ si , s 2 ) was calculated only recently [l5[ for the case of massless 
quarks. 2 The explicit expressions can be found in j 1 51 ] - 

The local-duality (LD) sum rules [H, H, correspond to the limit M — > 00. A remarkable feature of this limit 
is the vanishing of the condensate contributions to the sum rules for the pion form factor and the decay constant. 
Assuming the duality region in the si-s 2 plane to be a square of side so, and denoting the duality interval in the sum 
rule for the decay constant by sq, we obtain to a s accuracy 



flF^Q 2 ) = - 

T 

f = 



dsi 



ds 2 



A(°)( Sl , S2 ,Q 2 ) + a s A«( Sl , S2 ,Q 2 ) 
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~0(a 2 ) 
0(a 2 s ). 



(2.4) 
(2.5) 



One should not be confused by the simplicity of these expressions: the complicated nonperturbative dynamics is now 
hidden in the effective continuum thresholds so and so. Let us emphasize that the LD sum rules are predictive only 
if one knows, or fixes according to some criteria, the effective continuum thresholds. 
Some comments on the dispersion representations for f^F v and f 2 are in order: 



1 Notice that such a procedure of comparing two- and three-point sum rules finds a natural physical explanation within the correspondence 
between sum rules and the constituent quark picture observed in |17|1 . 

2 Another case in which the radiative corrections to the double spectral density of the three-point function have been calculated is the 
case of one massless and one infinitely heavy quark [HI , 
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1. Whereas the single dispersion representation for the decay constant f 2 is well-defined, the double dispersion 
representation for f^F^, even in the LD limit, has at least two essential ambiguities: 

a. The shape of the duality region in the S1-S2 plane: the simplest choice is a square, but any other region symmetric 
under si <-> s 2 may be also possible. 

b. Nothing forbids the upper boundary of the duality region so from being Q 2 -dependent, and additional assump- 
tions to fix so(Q 2 ) are necessary. Arguments in favour of choosing the parameters in two- and three-point sum rules 
constant and equal to each other were given in Let us see what happens if we choose the same constant value 
So = sq in the sum rules (|2.4p and (|2.5p . and substitute the sum rule (|2.5[) instead of f 2 into 



TT 2 J 


pso 

ds\ j 


pso 

' ds 2 




A (0) (s 1 , S2 ,Q 2 )+a s A( 1 )( Sl , S2 ,Q 2 )" 




pso 

' ds 




p(°)(s) + a.p«(a)" 





FA(f< - . ^ -• s^^^-. (2.6) 

The LD form factor given by (|2.6p has the following interesting properties [2(| : 

(i) It satisfies the normalization condition F^{Q 2 = 0) = 1 due to the vector Ward identity which relates the 
spectral density of the self-energy diagram and the double spectral density of the triangle diagram at zero momentum 
transfer: 

Jim A^(s u s 2> Q 2 )=npU( Sl )8( Sl -s 2 ), /»( s ) = -L p W( a ) = -L . ( 2 . 7 ) 

Q 2 ^0 47T 47T Z 

Clearly, for consistency one should then take into account the radiative corrections to the same order in two- and 
three-point correlators. 

(ii) Making use of the explicit expressions for AW ( one obtains 

^A<»>( S1>8S ,^) . te>, (2. 8 ) 



Substituting these expressions into (|2.6[) . one finds at large Q 



12. 



FAQ 2 ) = + + O (asft/Q*) + O (al) . (2.10) 

We find quite remarkable that the exact leads to the correct pQCD (up to the running of a s ) large-Q 2 asymptotics 
of the pion form factor obtained from the LD sum rule (|2.6p . Let us explain this important point: Whereas, e.g., 
the normalization of the pion form factor (|2.6p at Q 2 — is the consequence of the Ward identity, we do not see 
any rigorous condition which would guarantee the correct large-Q 2 behaviour when using the same value of the pion 
duality intervals in two- and three-point correlators. We find this to be a strong argument in favour of the universality 
of the pion duality interval. 

(iii) The O(l) contribution, shown in Fig.QJi, was calculated in [||: 

The explicit expression for the 0(a s ) contribution 

h(Q 2 ,s )= fds! [° ds 2 A^( Sl ,s 2 ,Q 2 ) (2.12) 
Jo Jo 

was obtained only recently [151 ]. Before that, it was proposed to use instead of the unknown integral the simplest 
Ansatz []]| 

which reproduces the value of the integral at Q 2 = 0, fixed by the Ward identity, and its asymptotic behaviour 
according to (|2.9p . Fig. [TJd compares the formula (|2.13p with the result of the exact calculation: as one can see, the 
proposed formula underestimates the exact jW by more than 20% in the broad range of practical relevance Q 2 =1-30 
GeV 2 . 
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Fig. 1: J^ji Io(Q 2 , so) (a) and Ii(Q 2 , so) (b) for so — An 2 f 2 . Solid (red) line: exact result, dashed (blue) line: Ansatz 

una). 



2. There are obvious problems with the application of this sum rule at small Q 2 < 1 GeV 2 : 



First, the OPE for the three-point correlator was obtained in the region where all three external variables \p±\, \p 21 



and Q 2 are large. Therefore, the sum rule cannot be directly applied at small Q 2 , although the expression (|2.6p leads 
to the correct normalization of the form factor. Additional contributions appear at small Q 2 [2lT | which prevent from 
giving unambiguous predictions in this region. Of course, allowing for a Q 2 -dependent value sq(Q 2 ) in the sum rule 
(12. 4p . we can formally extend the formula also to lower Q 2 and apply it starting from Q 2 = 0, but as we have noticed 
above, in this case the sum rule loses its predictivity. A technical indication that the LD sum rule (|2.6p cannot be 
applied at very small Q 2 is the presence of the terms ~ y/Q 2 [see (|2.11|) ] leading to an infinite value of the pion radius. 

Second, the spectral density A(si, S2, Q 2 ) contains the terms O(l) and 0(a s ), whereas higher powers are unknown. 
Since the coupling constant a s is not small in the soft region, our spectral density is not sufficient for application to 
the form factor at Q 2 < 1 GeV 2 . 

3. In order to apply the obtained formulas for large Q 2 , higher-order radiative corrections, leading to the running of 
a Sl should be taken into account. Such an accuracy is beyond our two- loop calculation; nevertheless, a self-consistent 
expression for the form factor applicable for all Q 2 > may be written as 

i fs (Q 2 ) Mo(Q 2 ) r -, 
F^Q 2 ) = ~7T~2 / ds i / ds * A^( Sl ,s 2: Q 2 ) + a s (Q* 2 )A ( - 1 \s 1 ,s 2 ,Q 2 ) , (2.14) 

JTT 71 ' JO JO 1 J 

where the scale Q* 2 in the argument of a s is related to Q 2 (see the discussion in (22|) and So(Q 2 ) satisfies the boundary 
conditions 3 

4tt 2 f 2 

s (Q 2 = 0)= ^ , So (Q 2 = co)=4 7 r 2 / 2 . (2.15) 

1 + a s (0)/7r 

If the effective threshold so(Q 2 ) satisfies these relations, the form factor is normalized to i^(0) = 1 and reproduces 
the pQCD asymptotic behaviour at Q 2 — » oo. 

In the following, we shall set the scale Q* 2 = Q 2 : in the region Q 2 > 1 GeV 2 , a s (Q 2 ) is a slowly varying function 
of Q 2 (Fig. [5^); therefore, the precise setting of the scale makes very little difference. We shall use the following 
appealing parametrization of s (Q 2 ), obviously satisfying (I2.15|) : 

4tt 2 f 2 

*o(Q 2 )= ^ ■ (2-16) 

Before turning to the numerical analysis, we would like to draw the reader's attention to the following observation: 
An essential feature of the form factor obtained from the three-point sum rule is the full cancellation of the double 
logarithmic terms. The proof of this general property of the color-neutral three-point Green functions in QCD can 



3 In [l9j it was argued that so is the relevant scale of a s in the LD sum rules for the decay constant and for the form factor at Q 2 = 0. 
For our discussion this subtlety is irrelavant so we somewhat symbolically use the notation a s (0). 
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be found in [3[ (see als o 12311 ); the cancellation of double logs was checked in explicit two-loop calculations for various 
quark currents in [Til . Il5j. In contrast to this result, the pion form factor obtained from the light-cone sum rule 
contains double log terms [HI]. This discrepancy requires a clarification. Presumably, higher- twist contributions, 
which were not taken into account in [l3|, but which are in general not suppressed compared to the lower-twist 
contributions [24| play a crucial role here. 



3. NUMERICAL RESULTS 



For numerical estimates, we make use of the three-loop running a s (Q 2 ) (Fig. [2^,). The corresponding sq(Q 2 ) is 
shown in Fig. [2b- Notice that it is a slowly varying function in the region Q 2 > 1 GeV 2 , where we apply the LD sum 
rule. 

The pion form factor is shown in Fig. [3^. The 0(1) and 0(a s ) terms, separately, are given in Fig. [3Jd. It should be 
noticed that the 0(1) term providing the l/Q 4 power correction at large Q 2 , dominates the form factor at low Q 2 , 
and still gives 50% at Q 2 = 20 GeV 2 . The 0(a s ) term gives more than 80% of the form factor only above Q 2 = 100 
GeV 2 . Such a pattern of the pion form factor behaviour has been conjectured many times in the literature; we now 
obtain this behaviour in an explicit calculation. The earlier analyses of the pion form factor in a broad range of 
momentum transfers [lU, [TH, [23, Hi| are consistent with the results reported in Fig. [3] within about 20% accuracy. 

Fig. H] presents the ratio of the O(l) and the 0(a s ) contributions to the pion form factor vs Q 2 for different models 
of the effective continuum thresholds. One can clearly sec that the ratio is mainly determined by the corresponding 
double spectral densities A^ ) and A' 1 ), whereas its sensitivity to the effective continuum threshold is rather weak. 
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Fig. 2: The perturbative a a (Q) (a) and the corresponding effective threshold so(Q 2 ) (b) given by (|2.15|l . Dashed lines show 
these quantities outside our working region. 
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Fig. 3: The pion form factor at Q 2 > 0.5 GeV 2 . Experimental data from [25[. Solid (red) line: the result of the calculation 
according to (|2.14|) . (a) Short-dashed (green) line: the form factor obtained with constant so = 0.65 GeV 2 ; long-dashed (blue) 
line: so = 0.6 GeV 2 . (b) Short-dashed (black) line: the O(l) contribution, long-dashed (blue) line: the 0(a s ) contribution. 
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Fig. 4: The ratio of the O(l) and 0(a 3 ) contributions to the pion form factor R(Q 2 ,so) = Io(Q 2 , so)/[a s (Q 2 )Ii(Q 2 , so)]. 
Solid (red) line: the result of the calculation with the effective continuum threshold so(Q 2 ) (|2.16p . short-dashed (green) line: 
so = 0.65 GeV 2 , long-dashed (blue) line: s = 0.6 GeV 2 . 

4. DISCUSSION AND CONCLUSIONS 

We have presented the analysis of the pion form factor in a broad range of spacelike momentum transfers making use 
of the local-duality sum rule. This is the first analysis which takes into account both the leading order O(l) contribution 
to the pion form factor and the recently calculated first-order 0(a s ) radiative correction. These ingredients are 
crucial for the possibility to consider the form factor in a broad range of Q 2 and to study the transition from the 
nonperturbative to the perturbative region. 

Let us summarize the essential ingredients, the uncertainties, and the lessons to be learnt from our analysis: 

• The double spectral density of the spectral representation for the form factor. We have good control over the 
spectral density — we have included the exact O(l) and 0(a s ) terms, and omit the (unknown) 0(a 2 ) terms, 
which are expected to contribute less than 10% at Q 2 > 1 GeV 2 . [The inclusion of the 0(a 2 ) terms in the 
spectral density would lead to a corresponding modification of the effective continuum threshold, with the net 
effect upon the form factor of only a few percent.] 

• The model for the effective continuum threshold: This very quantity determines to a great extent the value of the 
form factor extracted from the sum rule. The possibility to fix this threshold is the weak point of the approaches 
based on sum rules, which limits their predictivity [27| . 

We use the same universal effective continuum threshold in two- and three-point sum rules. This allows us to 
relate the value of the threshold to the pion decay constant, known experimentally. We therefore have no free 
numerical parameters in our analysis. 

There are at least two arguments in favour of our choice of sq(Q 2 ): 

First, we have demonstrated that it leads to the correct asymptotic behaviour of the pion form factor at Q 2 —> oo. 

Second, we expect our approach to work better with the increase of Q 2 . We have seen that it works very 
well already at relatively small Q 2 = 1-4 GeV 2 (recall that we have no numerical parameters to be tuned to 
reproduce the data). Therefore, we believe that for all Q 2 > 1 GeV 2 we give reasonable predictions. However, 
we cannot control the accuracy of our predictions for the form factor and cannot provide any error estimates. 

• We can, however, control much better the relative weights of the O(l) and 0(a s ) contributions to the form 
factor: their ratio is practically independent of the model for the continuum threshold and is determined to great 
extent by the corresponding 0(1) and 0(a s ) double spectral densities. Here, our results convincingly show that 
the 0(a s ) contribution to the pion form factor stays at a level below 50% at Q 2 < 20 GeV 2 and demonstrate in 
a largely model-independent way that the pion form factor is mainly of nonperturbative origin up to very high 
Q 2 . Thus, our results definitely speak against the pQCD approach to form factors at intermediate Q 2 , referred 
to as Scenario A in the Introduction, and confirm Scenario B. Although obtained without any reference to the 
shape of the pion DA, our results indirectly restrict the pion DA at low values of the renormalization scale: For 
instance, convex DAs of the type of [l2j], close to the asymptotic one, provide the form factor compatible with 
the results reported here. Also a broader class of the DAs, such as, e.g., a double- humped DA with a suppressed 
end-point region of [12] seems to lead to the pion form factor in agreement with our results. For a conclusive 
clarification of this point, the analysis of both the O(l) and 0(a s ) contributions corresponding to this DA is 
necessary. 
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Finally, let us notice that the local-duality version, formulated and developed by Radyushkin and co-workers, in many 
cases has definite advantages compared to other versions of QCD sum rules: For instance, the standard three-point 
sum rules cannot go to large Q 2 because of polynomial terms, the results from light-cone sum rules depend on the 
light-cone distribution amplitudes. Of course, as we have already mentioned above, the numerical results for the form 
factor from the local-duality sum rules depend crucially on the model of the effective continuum threshold used for 
the calculations, but this shortcoming is shared by all versions of QCD sum rules [27j. In addition to this uncertainty, 
other versions of sum rules have uncertainties related to parameters not precisely known, such as the condensates 
and the distribution amplitudes. We therefore believe to provide the most complete analysis of the pion form factor 
available for the time being. 
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